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We propose a "channelization" architecture to achieve wide-band electromagnetically induced transparency 
(EIT) and ultra-slow light propagation in atomic 8 'Rb vapors. EIT and slow light are achieved by shining a 
strong, resonant "pump" laser on the atomic medium, which allows slow and unattenuated propagation of a 
weaker "signal" beam, but only when a two-photon resonance condition is satisfied. Our wideband architecture 
is accomplished by dispersing a wideband signal spatially, transverse to the propagation direction, prior to 
entering the atomic cell. When particular Zeeman sub-levels are used in the EIT system, then one can introduce a 
magnetic field with a linear gradient such that the two-photon resonance condition is satisfied for each individual 
frequency component. Because slow light is a group velocity effect, utilizing differential phase shifts across the 
spectrum of a light pulse, one must then introduce a slight mismatch from perfect resonance to induce a delay. 
We present a model which accounts for diffusion of the atoms in the varying magnetic field as well as interaction 
with levels outside the ideal three-level system on which EIT is based. We find the maximum delay-bandwidth 
product decreases with bandwidth, and that delay-bandwidth product ~ 1 should be achievable with bandwidth 
~ 50 MHz (~ 5 ns delay). This is a large improvement over the ~ 1 MHz bandwidths in conventional slow 
light systems and could be of use in signal processing applications. 

I. INTRODUCTION 

Electromagnetically induced transparency (EIT) 1 1 , 2], in which a "pump" field of laser light can allow a weaker "signal" 
field to propagate through an otherwise opaque atomic gas, has been inspiring a number of applications based on the underlying 
coherent interaction of laser light with atomic media. These include nonlinear optics at low light levels 01 an d ultra-sensitive 
magnetic field measurements |4|] . 

Of particular interest has been the recent observation of ultra-slow light (USL) |5|, |6|] in atomic gases, at group velocities on 
the order of 10 m/s, due to a steep linear dispersion in the index of refraction associated with the narrow EIT feature. This could 
allow for controllable true-time delay devices for classical light pulses, with applications in fiber-optic telecommunications 0] 
and radar signal processing |8]. Later extensions of the technique to stored light 0, [Toll (for several milliseconds) has also 
raised the possibility of quantum memory devices II ill . While the narrow frequency feature (below the natural linewidth and 
Doppler width) of EIT is one of its attractive features for precision applications |4], this has drawbacks in regards to delay and 
storage applications. Optical communications and radar processing typically desire ~ 20 GHz bandwidth. Similarly, single 
photon sources and other tools of potential quantum information technologies may emit photons over a broad band. In USL 
experiments to date, the width of EIT transparency window is much narrower. 

EIT and USL work best when the atom can be well described with a A energy level structure. The signal field is near- 
resonant with a stable state (which we label |1)) and a radiatively decaying excited state (|3)). The pump field is resonant with 
another stable state |2) and the common excited level |3). We consider two energy level schemes in 87 Rb, shown in Figfja). 
The schemes are labeled "A" (dashed, blue arrows) and "B" (solid, red arrows). The transparency and slow, distortion-free 
propagation of the signal pulse that we desire occur only when the frequency difference of the two lasers ui s — ui p matches the 
energy level difference between levels |1) and \2) to within the narrow EIT width. Frequency components of the signal outside 
this width are strongly absorbed and distorted. This width is directly proportional to the pump power and practical limits on the 
pump power (~ 10 mW/cm 2 ) limit it to ~ 1 MHz. 

The need for a wide-band, controllable true-time delay device inspired us to here consider theoretically a "channelization" 
geometry utilizing Zeeman shifts in the atoms (see Fig^b)). There are a variety of techniques available to spatially separate 
various frequency components of broadband light. Assume a broadband signal pulse (represented as a dashed, yellow arrow) 
propagating along the longitudinal (z) dimension is split in the transverse (x) direction, with a transverse displacement propor- 
tional to the detuning from some chosen central frequency. In our channelization geometry, this dispersed signal then enters 
the 87 Rb cell, which is illuminated by a co-propagating, monochromatic pump field. When particular Zeeman sublevels are 
used, as in Figflja), the two photon resonance condition required for slow light will be a strong function of magnetic field along 
the quantization axis z due Zeeman shifts of the levels. Thus, applying a longitudinal magnetic field (dotted, green arrows) 
with a linear gradient along x will allow us to achieve the conditions for EIT and USL for each frequency simultaneously. The 
components can then be recombined after passing through the cell, resulting in a true time delay device for a broadband signal. 

USL is a group (rather than phase) velocity effect, meaning it works by applying differential phase shifts to each frequency 
component in the signal pulse. Maintaining perfect two-photon resonance everywhere would result in no differential phase shift 
and thus, no delay. Therefore, one should choose the magnetic field gradient such that there is a small, varying detuning across 
the cell. By choosing this mismatch to be small enough that all components are within the EIT resonance one can obtain the 
linear frequency dispersion necessary for slow light and still maintain transparency. There is a direct trade-off such that the 
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FIG. 1: Schematic of wideband atomic slow light system, (a) The hyperfine structure of the Di line in Rb. The bare energy levels are 
shown as dotted lines. The excited state splitting is (27r)817 MHz, the ground state splitting is (2-7r)6.834 GHz, and the D\ optical transition is 
at A = 795 nm tlUl . A magnetic field along the quantization axis (z) splits the individual Zeeman sub-levels |i) (solid black lines) according 
to = fJsgp' *rrSp B z shown schematically in the figure for a positive magnetic field. We consider two EIT A schemes (Scheme A as 
dashed arrows(blue online), Scheme B as solid arrows (red online)). The thicker arrows indicate the pump field, with frequency ui p chosen to 
be resonant with the 2) <-> |3) transition with no magnetic field. The thinner arrows refer to the signal field, at some arbitrary frequency u s 
within our desired bandwidth. The magnetic field is chosen such so the system is in two-photon resonance. In each case the atoms are initially 
purely in |1). In Scheme A, the fields also couple to |4) = \F' = 2, mf = —1), causing a decoherence mechanism which we investigate, (b) 
The channelization architecture. A signal field is split, such that the x position of each frequency component is proportional to the frequency, 
and then input into a cell illuminated by a monochromatic pump field. A magnetic field with a linear gradient (dotted arrows, green online) is 
then applied across the cell in such a way that two photon resonance is nearly maintained everywhere. 

delay -bandwidth product cannot be increased with this method. But for many applications, such as delays for signal processing, 
delays in conventional slow light systems are much longer than necessary, while the bandwidth is much too small. Our method 
allows one to circumvent this problem. 

To be of practical interest, such a system would ideally provide uniform transmission and delay over the entire bandwidth. 
We will see that Doppler broadening due to thermal motion of a room temperature vapor is actually beneficial to our scheme in 
this respect. Two-photon resonance is maintained for atoms of all velocities because the two fields are co-propagating, while the 
one-photon detuning strongly depends on the velocity of each atom. Once one averages over the Doppler profile one finds that 
delay and transmission insensitive to the one-photon detuning when it is well within the Doppler width (^300 MHz in room 
temperature rubidium) lH2ll . This means that the one-photon detuning resulting from our Zeeman shifts (see Fig.^a)) does not 
effect our delays and transmissions. 

There are several issues to consider to optimize the system we propose here. First, the transmission (and maximum delay) are 
limited by decay of the coherence between the two ground states. This decay time is often governed by the time it takes the atoms 
to leave the interaction region with the pump (due to thermal motion) |6]. This decoherence mechanism can be mitigated by the 
addition of a buffer gas (such as helium) I6. ll3lll4ll which significantly reduces the diffusion of the 87 Rb atoms via collisions. 
Second, decoherence can also occur from transitions to levels outside the A system of interest (for example, the level |4) in 
Scheme A in Fig.^a)). The presence of a high pressure of buffer gas can actually worsen this problem by pressure broadening 
these unwanted levels. This means there is often an optimal intermediate pressure which balances these two considerations. A 
third consideration, unique to the channelization architecture, is the presence of a high magnetic field gradient. Because the two 
photon resonance condition is only satisfied over a small range of magnetic fields, and therefore a small transverse spatial region, 
the thermal diffusion of atoms will put them out of two-photon resonance and potentially cause a severe absorption of the signal. 
This will mean higher buffer gas pressures may be desirable for our proposal than in conventional slow light. 

In this paper, we present calculations of the transparency, delays and bandwidths for both Schemes A and B in 87 Rb atoms 
with a helium buffer gas. We use a model based on the linear response of the signal in the medium, taking into account the 
diffusion of the atoms, the linear gradient of the magnetic field, pressure broadening, Doppler broadening, as well as interaction 
with additional levels in the hyperfine structure. The calculation uses a semi-classical model of the evolution of the ensemble 
average of the atomic density matrix interacting with classical light fields, and a classical treatment of the light field propagation 
based on the slowly varying envelope (SVE) version of Maxwell's equations 0]. 

In the first part of the paper we consider the case of a spatially uniform (but arbitrary) magnetic field and present a systematic 
analysis and optimization of the transmission and delays (with regards to buffer gas pressure and 87 Rb density, etc.). In the 
second part, we then introduce a model to account for a magnetic field with a steep linear gradient. We then discuss how to 
best choose the slight mismatch in the two-photon resonance to maximize the delay -bandwidth product. We find that Scheme B, 
provides better overall performance, with delays ~5 ns over a bandwidth ~50 MHz. We find that the delay-bandwidth product 
decreases with bandwidth, reaching unity at around 50 MHz. However, it is more difficult to prepare the initial state for this 
scheme (all atoms in \F = 1, ttif = +1)) than it is for Scheme A (all atoms in \F = 2, mp = —2), which is easily initialized 
via optical pumping from the pump). While Scheme A provides worse overall performance, it should still provide a suitable 
system to improve the bandwidth over conventional systems. 
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II. UNIFORM MAGNETIC FIELD CASE 

We begin by analyzing "conventional" slow light ij Uloll . with a homogenous magnetic field and attempt to understand the 
delays, bandwidths and delay -bandwidth products achievable in both schemes A and B. Our model accounts for effects of an 
extra level |4) and the buffer gas and can also explore the dependence on the magnetic field. 



A. Model 



We model the signal and pump light fields classically and represent them with Rabi frequencies £l s ( p ) = — er 13 ( 2 3) • 
ts(p)E s (p)/Ti, where are dipole matrix elements, e s ( p ) are unit polarization vectors, and -E s (p) are tne slowly varying en- 
velopes of the electric fields E s r p \ (e 1 ^^) 2 "'^?)*) + c.c)/2. Here the wavenumbers k s ,k p w 2n/A, with A = 795 nm the 
wavelength, are taken to be equal. The field polarizations e s ( p ) are chosen to match the transitions (er+, er_ in Scheme A and 
<7+, cr + in Scheme B). Meanwhile the 87 Rb atoms are represented by a 4 x 4 density matrix p° (representing each internal level 
|i) under consideration). The diagonal elements represent the fractional populations in each state while the off-diagonal terms 
represent the coherences between levels (induced by the coherent lasers). At the microscopic level, each atom has a density 
matrix but we course grain average over spatial regions large compared to the inter-atomic spacing, meaning the total density 
matrix can be approximated by Np°, where N is the atomic density and p° is normalized to unity ^ p° = 1. 

Our model for the atomic density matrix evolution has two parts: (1) a coherent, Hamiltonian part, which includes coupling 
with the light fields, Zeeman shifts, and pressure shifts; and (2) an incoherent part, which includes spontaneous decay, pressure 
broadening, and diffusion out of the pump field interaction region (or the cell walls if the pump illuminates the entire cell). The 
evolution equations for the density matrix elements are given by: 
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The parameters /3i4(24) = e s ( p ) • i"i4(24) /^s(p) ■ r i3(23) characterize the difference in coupling (both sign and amplitude) for 
the unwanted transition to |4). They are given by ratios of Clebsch-Gordon coefficients and are /3i4 = 1/V3, /?24 = — \/3 for 
Scheme A 1 15] and vanish for Scheme B, where there is no coupling to additional levels. We have made the rotating wave 
approximation to eliminate counter-rotating terms and coupling to levels detuned by the ground-state hyperfine frequencies (~ 
6.8 GHz), but kept terms detuned only by the excited-state hyperfine detuning (~ 800 MHz). The bare level frequencies ujj on the 



diagonal terms are shifted by linear Zeeman shifts A^ = iisgp'rUp' B z where the Bohr magneton is p,s = (27r)1.4 MHz/G, 
the Lande g-Factors Qp are given in Fig.^a), and B z is the magnetic field 1 15]. They are also shifted by buffer gas pressure 
shifts Sp p, where p is the pressure. For 87 Rb with a helium buffer gas, these are taken to be S p e ^ = 
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(the ground state pressure shifts are much smaller and negligible for our 



for all excited manifold (F' 
parameters). 

The incoherent evolution is governed by a super-operator C. Writing out only the non-zero terms 
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C\l = -7diff, = ~( T r/ 2 + B vV), £\\ = ~(T r + 2B pP ) (3) 

(plus the corresponding complex conjugate terms for the off-diagonal elements). These terms are generally non-Hermitian as 
they involve lossy terms and incoherent transitions due to spontaneous emission. The first line represents feeding the ground 
states |1), |2) via spontaneous emission from |3), |4). The rate of emission from these levels is T r = (27r)6 MHz but this can 
branch into both levels |1), |2) as well as other levels outside of the system of interest. The branching ratios to various states 
are given by the oscillator strengths /y which are proportional to the square of the Clebsch-Gordon coefficients. For Scheme 
A, /13 = 1/2, /23 = 1/12 and fj4 = \(3ji\ 2 fj3, while for Scheme B, /13 = 1/2, /a 3 = 1/6. The second line represents the 
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population loss of the excited states from spontaneous emission. In the third line, we have put in the dephasing of the ground 
state coherence due to the diffusion of the gas out of the interaction region (taken to be of width u> int and height /i; nt ) I14ll : 

2 9 D a 1 
7diff = -2.405 2 f— — — — (4) 

O Wi a t h hlt 1 + 6 . 8 < m f p / V Wint hint 

where the diffusion constant for 87 Rb in a helium buffer gas of pressure p and temperature T is 

D g = (410 cm 2 /s)(Torr/p) ^T/273 K O. The mean free path is Z m f p = 3D g /vth and the thermal velocity v t h = 
jm. The factor (2/3) in front is to account for the fact that length (along z) is generally much longer than w; nt , h- lnt 
and diffusion in this dimension does not effect the coherence. There is additionally a depolarizing cross-section from 87 Rb-He 
collisions, which can dephase the ground states, but this is much smaller effect than other decoherence mechanisms for our 
parameters of interest and is neglected. The last two lines of Eq. (0 contain dephasings of coherences from radiation, at IY/2, 
and pressure broadening, with B p = (2tt)5 MHz/Toit II 8tll9ll . 

When solving Eq. Q, we consider the usual weak signal regime, and drop all terms higher than linear order in Vl s . This is 
valid when |f2 s | <C |fip| and when multiple scattering of spontaneously emitted photons can be ignored (which is usually the 
case in EIT since spontaneous emission is suppressed). Then we can take p^i — > 1 ar, d P221 /°33> /°44j P231 P241 P34 0- We 
are left with three non-trivial equations for the evolution of p^, p 31 , p 41 . We furthermore make transformations to eliminate the 
time-dependent terms in Eq. (|2}: psiai — P31 41 e~ l ( fcsZ_Ws *', p 2 i = ( o° 1 e~ l (( fcs ~ fe p) z ~( Ws ~ w J>)*). For convenience we define a 
vector p = [P21, P31, P4i] T - The evolution Eq. Q can then be written as: 



p = Mp + S; 

fi(A' s -A' p )- 7diS -in; 
m = -%n p i(K + 5 D ) 

S EE ( ,-^ s ,-^ 14 ft s ) T 

where we have defined the shifted detunings, A' s = A s - A^ - S p c) p + A^\ A' p = A p - - S p c) p + A^ , and 

A 43 = ujf^ + A% — lu^ — A% . The bare detunings are A s ( p ) = u> s ( p ) — (ll>3 — w i(2))- Throughout we choose pump to 
resonant with the bare resonance A p — 0, while the signal A s varies. We have also introduced a Doppler shift 8d — (27r)u z /A, 
where v z is the velocity of a particular atom along the light propagation direction z. The total dephasing rates of the optical 
transitions are 7 e = IY/2 + B p p. 

When studying the light field propagation, it will be easiest to work in Fourier space and so we transform from p{t) — > p(8) 
and n s (t) — > ft s (5). Equation Q is linear in time-dependent quantities and finding the solution of its Fourier transform is 
equivalent to solving for its steady state in the time domain but replacing A s — > A s = A s + <5 where <5 represents the deviation 
of a particular Fourier component of the signal field ft s from the central probe frequency lo p (due to time dependence). This 
solution is: 
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p = A4- 1 S, (6) 

where M^ 1 is simply Ai^ 1 after making the replacement A s — * A s and in S we replace il s — * Ct s . 

Finally, to obtain the response of the entire medium we integrate over the thermal profile of velocities v z 01 : 
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p {D) (A s )= I dS D p(A s ,5 D )exp[ ~-^r), (7) 



where the Doppler width is Ad = -\/2/3(27r)w t ^/A (the square root term is a geometrical factor). 

Turning now to the light propagation, in the linear signal regime the pump field Rabi frequency il p is constant in space and 
time. The signal field ft s propagates according to the S VE Maxwell equation, with the polarization written in terms of the atomic 
density matrix. Once we Fourier transform the Maxwell equation, it can be written: 

—n s = l Nf 13 a X {D) (& s )n s + i-ri s 

oz 2 c 



wVro X W(A S ) = +^ 14 p(f) (8, 
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FIG. 2: EIT transparency window and steep index of refraction, (a-b) Frequency dependence of the real (a) and imaginary (b) parts of the 
susceptibility X m Schemes A (thin curves, red online) and B (thick curves, blue online) with B z = 0. Buffer gas pressures are p=3 Torr 
(solid curves) and 30 Torr (dotted curves). The scale is too large to see the EIT resonance features, which is located at the vertical dotted lines, 
(c-d) The same curves shown on a smaller scale, at the EIT resonance. The dotted lines indicate exact two-photon resonance A' p = A' c . 



and a = 3 A 2 /(2tt) is the resonant cross section for unity oscillator strength. The last term accounts for the propagation at c in 
free space. If the fastest time scale of interest in the problem is much slower than total length of propagation (~ cm) divided by 
c then this term can be neglected. Note that since /O31, /?4i oc Cl s , the susceptibility is independent of Cl s . Equation l|8) can be 
trivially solved to give how a particular frequency component, with an input amplitude f2g (A s ), will propagate through a cell 
of length Z co n in the medium f2s° ut ^(A s ) = ei Dx(D) ^^fli in \A s ), where the "optical density" is defined as D = fi3Nal ce \\ 
(we have ignored the speed of light propagation term for simplicity). From this one clearly sees how the imaginary part of x^ 
is proportional to the absorption cross-section while the real part is gives phase shifts and thus determines the index of refraction. 

In a non-Doppler broadened, non-pressure broadened medium in the absence of a pump field (il p — > 0), and with no level |4) 
(/3i4 = 0), we recover the usual Lorentzian susceptibility profile. The absorption cross-section is peaked at the atomic resonance 
A' = and has a width r r while the real part exhibits anomalous dispersion at the resonance. The dimensionless susceptibility 
X- D ^ in Eq. (|8} has been defined such that the peak absorption is x^ (Ag = 0) = i, leading to exponential attenuation of the 
signal intensity I s oc e~ D . Pressure broadening and Doppler broadening act to reduce this resonant cross section while widening 
the feature. 



B. Results 



The presence of a sufficiently strong pump field Q, s will then introduce a sharp EIT feature at frequencies very near the 
two-photon resonance, in the middle of this broad absorption resonance. At the EIT resonance, we have reduced absorption 
(transparency) and a linear slope in the index of refraction with normal (positive) dispersion (leading to the slow group velocity). 
We now proceed with some calculations of these EIT features. The primary motivation of the present section is to learn the role 
played by the buffer gas pressure and to quantitatively learn the dependence on magnetic field. These results will later guide the 
optimization of the channelization architecture with the linear magnetic field gradient. 

Throughout this section consider a case with a 1 mW pump laser with an area tOmt/Mnt = (2 mm) 2 . Due to differing oscillator 
strengths ^23 this results in Rabi frequencies of Cl p = (27r)8.45 MHz for Scheme A and il p = (27r)12.3 MHz for Scheme B. At 
a temperature of T — 333 K (60 degrees Celcius), there is a 87 Rb density of N = 2.5 x 10 11 cm" 3 [15] and a Doppler width 
A D = (2tt)320MHz. 

Considering first a homogenous B z = magnetic field, we plot in Fig|2]the real and imaginary parts of the Doppler averaged 
susceptibility x^ (A s ), x\°^ (A s ), according to Eq. (|8}, both on a large (a,b) and small (c,d) scale for both Schemes A (thinner, 
blue curves) and Scheme B (thicker, red curves). From Fig|21a,b) we see that, away from the narrow EIT feature, the susceptibil- 
ity retains the usual Lorentzian susceptibility feature one would expect in a two-level medium, with the width ~ Ad and height 
~ T t /2Ad- The biggest apparent difference between the two schemes is that in Scheme A there is an extra resonance, due to 
level |4), near the hyperfine splitting (~ 817 MHz). Note that we are in a regime where the Doppler broadening is comparable 
to this splitting, causing the two resonances to slightly overlap. The solid curves are for the relatively small buffer gas pressure 
p — 3 Torr, while the dotted curves are for p — 30 Torr. At this higher pressure, the pressure broadening becomes important as 
7 e = (2tt) 253 MHz becomes comparable to Ar>. 

Examining the smaller scale plots Fig.|2jc-d), we see the pump field introduces a sharp feature at the two-photon resonance 
(A s = 0). In the imaginary part we see a narrow notch in the absorption cross-section (the transparency window) while in the 
real part we see steep linear dispersion, the origin of the slow group velocity. To describe the degree of transparency, we define 

the parameter -Reit as the ratio of the minimum of Xi (A s ) in the resonance to its value there with no pump field present. We 
note a very large difference in this parameter for the two schemes, with only Kbit ~ 1/3 for Scheme A. As we will discuss 
below, coupling to level |4) is the reason for the lack of good transparency in this case. We also note that the transparency 
window is slightly shifted from exact two photon resonance (by approximately —0.1 MHz), an effect of stark shifts from |4). In 
Scheme B, the transparency is quite good i?EiT ~ 0.001. 
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FIG. 3: Pulse propagation characteristics (a) For the case p = 2 Torr, B z = in Scheme A we calculate A = 0.005, S = 24 ns, and 
W = (2-7r)0.72 MHz. Here we plot the resulting delay td = SD/2 (solid curve, blue online), loss L = (l — e~ AD )) (dashed, red) and inverse 
bandwidth = \f~D/W (dotted, green) versus optical density D. (b) Same for Scheme B, withp = 15 Torr, for which A = 3.5 x 10 -5 , 
S = 13.3 ns, and W = (2vr)2.31 MHz. 



To translate these x (A s ) curves into performance for delay devices, we examine the shape of these curves at the EIT 
resonance. To a good approximation, the real part here can be written as some central value plus a linear part, while the 
imaginary part is some minimum value plus a parabolic shape: 



X W (A s )^o + S(A s -Ao) + »A + » l s ^ °> (9 ) 

where Ao is defined to be the minimum point of Xi an d the parameters 0o , A absorption S (slope) , W (width) are ob- 
tained by numerically evaluating x^ (A s ) and its derivatives at Ao. In these terms our transparency parameter is -Reit — 
A/Im{x^(A )\n p ^o}- 

Now suppose we input a signal pulse a central frequency A s = A and a 1/e intensity half-width r s (giving it frequency 
components in the range A s = A ± t^ 1 ). To calculate what the pulse looks like after propagation through the medium, 
one takes the Fourier transform of the input pulse, calculates the propagation of each Fourier component according to Eq. (jSJl, 
then inverts the Fourier transformation back to the time domain. If r s sufficiently long that Eq. l|9} is valid for all frequency 
components of the pulse (which is usually true when t" 1 <§; W), then one can analytically perform the inverse transformation. 
The result is a pulse which is delayed in time by td — DS/2 and attenuated by a factor e~ AD . This Fourier analysis of 
the propagation reveals clearly how the delay comes about by a differential phase shift of the different frequency components 
SA' s D/2. One immediately sees there is a trade-off between transparency and delay. In addition, frequency components 
slightly off the resonance will be preferentially absorbed. This leads to a time-broadening of the pulse by r — > \/t| + (3~ 2 , 
where (3 — W/\/~D, and a reduction in the peak intensity by t 2 /(t 2 + (3~ 2 ). In this sense, (3 can be interpreted as the bandwidth 
of the system, as 8 ~ t^ 1 <C is required to prevent attenuation and distortion of the pulse. 

An obvious figure of merit is the delay-bandwidth product, important for many applications such as quantum memory devices 
or optical buffers. This parameter indicates the degree of "pulse separation" one can achieve (i.e. the number of pulse widths 
one can delay). Incidentally, this parameter being unity corresponds to the point at which the differential phase shifts across 
the frequency spectrum of the pulse is ~ (2tt). However, the absolute value of the bandwidth and maximum delay are also of 
importance, depending on the application. All these parameters depend strongly on the optical density D, which can generally 
be adjusted since the 87 Rb density N is a strong function of temperature 11511 . At N = 2.5 x 10 11 cm -3 for a i ce n = 1 cm cell, 
D = 370. For a given desired transmission, the parameter A tells us the optical density through which a pulse can successfully 
propagate. For a transmission 1/e we need D < D max = A^ 1 . This, in turn, determines the maximum achievable delay 

T^max) = S/2A. The delay -bandwidth product is td(3 = SW\f~D /2. This product increases with D and the best possible 
delay -bandwidth (using D max ) is SWj2\J~A. 

For Scheme A we find the performance is better for lower pressures. At p = 2 Torr and other parameters as in Fig. [2] we 
calculate A = 0.005, S = 24 ns and W — (27r)0.72 MHz, giving a maximum delay-bandwidth product of 0.76. In Fig. 13 a), 
we plot resulting pulse propagation parameters versus D. We see that for D ps 50 one obtains a to = 0.5 /is delay with only 
20% loss. But the minimum pulse width here is = 1.5 /is > r™ ax One does not get into the "pulse separation" regime until 
at about D = 350, where td = fi^ 1 = 4.1 /is where the attenuation is 82% and the bandwidth is (3 = (27r)38 kHz. 

By contrast, we can do much better in Scheme B. We find this system is better at higher pressures and for p = 15 Torr case 
we calculate A = 3.5 x 10~ 5 , S = 13.3 ns, and W — (27r)2.31 MHz, giving /Jr^ 11 ^ = 31. The pulse propagation parameters 
are plotted in Fig.|3lb) and we see the pulse separation regime begins at D w 100, at which point td ~ 0.7 /is and the loss is 
only L w 0.003. The bandwidth here is (3 = (27r)230 kHz. Note that despite the better performance, one would still have to go 
to very small D to get the bandwidth significantly more than 1 MHz, and here the delay -bandwidth product is quite small. 

We now attempt to get some intuitive understanding of what determines A, S, and W in the two schemes and study the 
dependence on buffer gas pressure and magnetic field. Scheme B, because of the absence of level |4) is significantly easier to 
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FIG. 4: EIT in a homogenous magnetic field. Plot of \ m the two schemes, again for p = 3 Torr. The solid curves are for B z — 
(identical to Fig. [2}, while the dotted curves are with a magnetic field B z = 75 G, which shifts the two-photon resonance (A' s = A' p ) to 
s ICB B z = (2tt)105MHz. 

understand. One can analytically obtain the solution Eq. and plug it directly into the susceptibility in Eq. (|8|i). The conditions 
necessary for EIT are that the ground state decoherence rate is small •y^is <C 7 e and the pump field intensity is sufficiently strong 
|f2 p | 2 3> 7diff7e- Assuming these inequalities and Taylor expanding in A' s about one- and two-photon resonance (A^ = A' p = 0) 
one obtains a susceptibility in the form of Eq. l|9} with the parameters cf>o = 0, A = 27diffr i ./|f2 p | 2 , S — 2r r /|f2 p | 2 , W — 
\Qp\ 2 / V^le^r- Thus we find that the absorption scales as the inverse of the pump intensity, the bandwidth scales directly with 
the intensity and the delay scales inversely with the intensity. There is a trade-off between delay and attenuation in choosing 
the pump intensity (just like the optical density). However, unlike with optical density, one cannot improve the delay-bandwidth 
product by changing the pump intensity. One can increase the bandwidth (and reduce the delay) but extremely high bandwidths 
require unreasonably high pump powers. We find that these analytic expressions are still valid after Doppler averaging. For the 
case we have just considered (Fig.|3lb)) these estimates yield A — 3.5 x 10~ 5 , S = 13.3 ns, W = (2-7r)2.34 MHz, in excellent 
agreement with the numerical results. 

At the point of transparency we find the ground state coherence is P21 ~ —Q s /tt p , which is known as a "dark state" Q]. The 
strong pump field acts to drive the system into this state. In this case the two terms driving absorption, on the |1) <-* |3) and 
1 2) 1 3) transitions, (see the first two terms on the third line of the Hamiltonian, Eq. (|2j), are equal and opposite, leading to a 
quantum interference which suppresses the absorption process. However, a non-zero detuning A' p or decoherence 7diff causes 
P21 to slowly evolve out of the dark state. Our analytic expression A = 27diffr r /|f2 p | 2 reflects the steady-state which occurs 
due to the balance of the preparation of the dark state by the pump field and the loss from it due to diffusion. To minimize our 
absorption, we clearly want to minimize 7diff, which can be done by increasing the buffer gas pressure (which decreases D g ) 
or increasing the effective area of interaction h lrA x u>; nt (see Eq. @). In practice, there is a trade-off between the interaction 
area and \Vl p \ 2 as one can increase the pump intensity by focusing the it more tightly. Thus, for a given pump power and buffer 
gas pressure, there is only a marginal dependence of A on the focusing area. Numerically we fine that there is a slight benefit in 
tighter focusing. Regardless one can increase the pressure to improve A to the desired level. The slope S is almost completely 
unaffected by p, while W decreases with pressure due to the factor ^/t7 in the denominator. 

For Scheme A, the situation is significantly altered by the presence of |4). The problem is that the dark state with respect to 
absorptions into |3) (the state for which the two absorption channels have equal and opposite amplitudes) is p2i = —Cl s /Ct p 
while the dark state with respect to |4) is P21 = —0i4^l s /024^p- Thus, unless oscillator strength ratios are equal (3\4 = 024 
there will be some absorption present for any value of ,021- In Scheme A, 0\4 = and 024 = ~ l\/3- 

This problem was studied in detail for the cold atom (non-Doppler broadened) case in 12011 . There it was found that, for 
atoms nearly resonant with |3), this effect leads to minimum absorption coefficient A g _ ros = (0\4 — /324) 2 7err/4(A 4 2 + 7 2 ) 
and an AC Stark shift of the EIT resonance by Aac = 024(024 — /3i4)|^p| 2 7e/4(A4 2 3 + 7 2 ). In the Doppler broadened case 
this problem is further complicated by the fact that, when the Doppler width Ac is comparable to the hyperfine splitting A 43 , 
there is a significant fraction of atoms which interact with both |3) and |4) with similar strength. The EIT interference is almost 
completely destroyed for these atoms. Numerically, we find that for Scheme A, this increases j4 ff-res by nearly an order of 
magnitude from the analytic estimate above. Increasing pump power increases the coupling to |4) in such a way that it exactly 
offsets any gain in the strength of the EIT resonance, and so we find Rbit quickly saturates to a value around 0.3 (as in Fig. [3c- 
d)), a constant determined by the relative values of Ad and A 43 . The interesting thing to note about y4 ff_ rcs is that it increases 
with buffer gas pressure (via j e ) because pressure broadening increases the relative role played by coupling to |4). Therefore, 
the diffusion problem favors higher p, while the interaction with |4) generally favors lower pressure and the lowest overall 
absorption is achieved by balancing these two considerations. We find that the optimal pressure is rather low (between 0.5 - 2 
Torr, depending on interaction area), but that the dependence is rather weak and so i?EiT ~ 0.3 is a good estimate over a broad 
range of pressures and interaction areas. The AC Stark shift is visible in the plots Fig.|5Jc-d) and agrees well with the above 
expression. 

Finally, in preparation for our channelization calculations, we must consider the degree to which our results depend on 
the value of the homogenous magnetic field B z . We still choose our pump field to be resonant with the bare |2) «-► |3) 
transition A p — oj p — (u>3 — 0J2) = 0. This system will then be in two-photon resonance for a probe photons with bare 



8 




B, (G) 



0.014C 
0.012 
0.01 

Ul o.oos 

3; 0.006" 
fn 0.004 
0.002 



B 2 (G) 




B z (G) 



FIG. 5: Dependence of EIT resonance on magnetic field. The parameters A (a) S (b) and W (c) , versus magnetic field B z . For Scheme A 
(thinner, blue), the pressure is p = 1.2 Torr, and for Scheme B (thicker, red) p = 30 Torr. 



detuning A s = A^' — A^' — Aac = SrcsB z — Aacs where s ros = VB{g^p — gpirip). In both Schemes A and B 
s ros = (27r)1.4 MHz/G. In Fig.|4]we plot x f° r both B z = (solid curves) (as in Fig.|2j and with a large magnetic field 
B z = 75 G (dotted curves). In Fig.[4la-b) we see that the magnetic field shifts the overall one-photon detuning so the Doppler 
broadened Lorentzian resonances are shifted. It also shifts the two-photon resonance (vertical lines) by a slightly different 
amount (so the EIT resonance with B z = 75 G is not exactly at peak of the Lorentzian resonance). However, this difference is 
still well within the Doppler width. Examination of the gp and rap shows that the signal detuning from |3), A^, (see Eq. (|5} and 
below) is only 1/6 of that of the two-photon resonance shift s les B z (this is true in both Schemes A and B, though the relative 
sign of the shift for the two schemes is opposite). Thus, even a shift of the two-photon resonance s 1QS B z ~ Ad will result 
in substantially smaller one-photon detuning. This is beneficial as the widths and strengths of the EIT resonances are effected 
when the one-photon detuning becomes comparable to the Doppler width 1 12]. Comparison of plots of the EIT resonance in 
Fig.0Jc-d) with the B z = case in Fig.|2jc-d) reveals they looks almost identical. 

To check numerically that the EIT resonance is indeed insensitive to B z over a wide range, we calculated the parameters A, S 
and W over the range B z =-200 G to +200 G. The results, presented in Fig.|5]bear out our expectation. In Scheme B there is no 
visible dependence on B z on the scale plotted. There is a weak dependence in the parameters in Scheme A, again primarily due 
to the effect of level |4) . The total range of resonance frequencies for this range of magnetic fields is (27r)560 MHz. 



III. EIT AND SLOW LIGHT WITH A MAGNETIC FIELD GRADIENT 



We are now prepared to add a final piece of our model account for the large gradients in the magnetic field and the transverse 
spatial dispersion of the signal in the proposed channelization geometry. We will then use this to characterize the performance 
of the delay device in each scheme, and explore quantitatively the maximum delays versus bandwidth and the optimal buffer gas 
pressures. 



A. Model 



Let us first consider how the susceptibility is affected by the diffusion of the gas in the presence of a magnetic field gradient. 
When a gas diffuses with some diffusion constant, it's density matrix evolves as p = ■ ■ ■ + D g V 2 p 1151 so we add this term to 
our existing evolution equation. Since the magnetic field gradient only exists in the x dimension, the diffusion in y and z has no 
effect. We then write Eq. (|5j: 



p = Mp+S + ^-^p. (10) 

With a magnetic field gradient our shifted detunings A^, A^ becomes spatially dependent and we take -j^p — * (^fjf") 2 ^rg^P- 
Suppose we apply a linear gradient so B z (x) = Sbx over the interaction region of width Wi n t- Then the two-photon resonance 
would vary by f3 c s = SsuimtSres over m e width of the interaction region, determining the effective bandwidth of our system 

Our susceptibility in the presence of diffusion will then be the steady state solution of Eq. (1101 . It is difficult to find this 
analytically in general, but this can be achieved with a perturbation approach under the assumption that the new diffusion term is 
small. In the perturbative approach, we take a zeroeth order solution to be solution without the diffusion term p^ — A? _1 S 
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FIG. 6: Schematic of detunings in wideband slow light scheme. The signal pulse frequencies are dispersed along x with a slope S'diap, 
indicated by the black solid arrow, which is slightly larger than the slope of the two-photon resonance Sbs ibs , indicated with the dashed arrow. 
The values chosen for the mismatch in slopes is such that the detuning local resonance A' p varies linearly with x and is everywhere within 
the local bandwidth /3. The top row diagrams this mismatch with respect to the local resonance at several locations, with the solid and dotted 
curves showing, respectively, the general shape of the c/ii^ 3 ' , x\ D ^ at each location, and the crosses indicating exact two-photon resonance. 



(see Eq. (|6}), and then plug this back into the full equation to obtain a correction from the diffusion term: 

(U) 

Because the two-photon detuning A^ — A^ is directly proportional to B z , the second derivative with respect to the B z basically 
corresponds to the curvature with frequency of at the EIT resonance, which Eq. (|9jl predicts to be 2/W 2 . In a sense, the 
diffusion term causes an averaging over some frequency width, which will partly wash out the EIT resonance, increasing our 
minimum absorption A. It also has a tendency to widen the feature, increasing W and decreasing our S. 

After obtaining our corrected density matrix p = p^ ' + p^ we once again Doppler average according to Eq. Q and then 
calculate the susceptibility. Unfortunately, the perturbative procedure is not valid near the wings of the Doppler profile, where 
the EIT feature becomes very narrow, even in cases where the perturbation is small near the center of the Doppler profile. In the 
limit that the diffusion correction becomes large p should smoothly return to it's value without the EIT feature p\n p ^o, but in 
Eq. (lilt the correction p~W can grow without bound. For this reason, we must use a slightly more complicated procedure, which 
in the limit of a small diffusion term reproduces Eq. (f^i and in the opposite limit reverts to p|o _, . This is accomplished by 
taking an average over a small range of magnetic fields: 

p(A s ,S D ,B) = J dS B ^ K -p (0 HA s ,S D ,B z +5 B )Exp^ - 



where A 2 B = 



4D g S% 


Ej V'; v ; 


3 


E, 'V, 



(12) 



and the and {vj} are, respectively, the eigenvalues and eigenvectors of M. 1 and the cij are the coefficients aj = Vj ■ 
(<9 2 p (0) /dBf). After calculating Eq. O we can then Doppler average with Eq. 0. 

Now to consider the spatial dispersion of the signal light, suppose that the signal frequency is varying linearly with x: A s = 
SdispX. If we chose Satsp = Sbs IC s then we would be in perfect two-photon resonance everywhere. However, this results in no 
differential phase shift across the spectrum of the pulse and thus, no delay. To achieve a delay, one must choose the these two 
slopes to be slightly mismatched, as diagrammed in Fig. [6] So long as the local detuning at the edges of the cell (or interaction 
area) \A' S \ = (Sdisp — <S , ss r es)(0.5wi n t) < /3 = W/y/D, the signal light is everywhere locally in the EIT regime. Then we 
recover a differential phase shift across the spectrum of the pulse, characterized by the effective slope: 



Pdisp JBSros) q 

OeS = ^ O. (13) 

»^disp 

If one chooses the maximum allowed mismatch £( max ) then ratio by which are delay decreases S c s / S is the inverse of the ratio 
by which our bandwidth increases /3 c g / j3. 

We note that the AC Stark shift in Scheme A must be accounted for to properly choose the mismatch in frequencies. However, 
the differential AC Stark shift is linear with frequency and so can be compensated for if needed. 
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FIG. 7: Effect of magnetic field gradient on the EIT resonance, (a-b) The real (a) and imaginary (b) parts of x f° r Scheme A with 
a pump power 5 mW and interaction area Wi n t — 2 cm, Zj n t = 0.5 mm, and pressure p = 25 Torr. The solid (blue online) curves show 
the resonance without any magnetic field gradient, while the black and gray dots are calculated with Sb = 2 G/mm and Sb = 8 G/mm, 
respectively, at the point where B z — 0. (c-d) The same calculation for Scheme B (red online). 
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FIG. 8: Dependence of EIT resonance on gradient, (a) The relative EIT absorption Reit = A/x (^o)\n p ->o versus magnetic field 
gradient in Scheme A (smaller, blue dots) and Scheme B (larger, red dots). We use a pump power 5 mW and interaction area i^mt — 2 cm, 
lint =0.5 mm and a pressure p — 10 Torr for Scheme A and p — 25 Torr Scheme B. (b) The slope S of the resonance for the same parameters, 
(c) The width of the resonance W. 



B. Results 

In Fig. we present examples of local susceptibilities at x = (B z — 0), calculated with the above procedure in both 
Schemes A and B and with various magnetic field gradients. The solid blue and red curves show the case with no gradient while 
the black and gray dots show the results with the gradients Sb = 2 G/mm and Sb = 8 G/mm, respectively. As expected the 
higher gradients wash out the EIT resonance, reducing both the transparency and the slope. This leads to a natural trade-off 
between obtaining a higher bandwidth /3 e g (with larger gradients) and better transparency (with lower gradients). As in the 
homogenous case, Scheme B offers better EIT for any given Sb, though it is worth noting that the transparency in Scheme A, 
because of the problem already present with level |4), is much less sensitive to the introduction of gradients. 

An analytic treatment calculating the correction Eq. using the non-Doppler broadened x obtained from Eq. (|6}, reveals 
that the expected absorption at the resonance (in the absence of other decoherence mechanisms from |4) and 7diff) is A gia d = 
64:(Dg/3)S B S^ cs T r ^/\flp\ 6 . We performed full numerical calculations of x , using Eq. d!2i . choosing values of h m t — 
0.5 mm and w lxx t — 2 cm and a pump power 5 mW, which gives tt p = (27r)12.0 MHz for Scheme A and fl p = (27r)16.9 MHz 
for Scheme B. In Fig.|8]we show the dependence of the absorption A, slope S, and width W versus Sb- The values chosen for 
the pressure correspond to optimal choices we discuss later. Our analytic expression above provides a reasonable estimate but 
is not quantitatively accurate. In particular, note the predicted quadratic dependence of the absorption A grac i ~ S B holds only 
for very small values of Sb then it becomes close to a linear dependence. As expected, the gradient contribution is dominant 
in Scheme B even for very small values of Sb, while in Scheme A the |4) contribution is dominant until about Sb = 4 G/mm. 
One also notes the gradient substantially impacts the slopes S. 

The performance as a delay device can be characterized again by the maximum delay with 1/e attenuation. For this analysis 
will choose the maximum optical density D max = 1/2 A and then demand that the maximum detuning at the cell edges, 
(Sdi S p - 5ijs r cs)(0.5u; in t) = <5 (max) = W/^2D max Q. This results in an effective slope S cS = SV2AW/l3 cS and a delay- 
bandwidth product (5' c ff-D max /2)/3 e g = WS/ V 2 A, virtually identical to conventional slow light. Examination of Fig.[8]shows 
that the increased absorption A and decreased slope S in fact decrease this product with bandwidth. 

In Figs. |5Ja-b), we show the effective dispersion curves for Scheme B the cases in Fig0c-d), with the mismatch chosen to 
be the maximum allowed difference as just described. These curves are just stretched in A s , with the stretch factor being the 
inverse of the ratio on the right-hand side of Eq. Jl 31 . The black and gray horizontal bars indicate the bandwidths /3 e g (the 
frequencies dispersed within the w- ln t = 2 cm area) in each case. In Fig. [Sic) we plot the maximum delay versus the bandwidth 
in both schemes. It appears delays ~ 5 ns are possible over a bandwidth ~ 50 MHz in Scheme B, whereas similar delays are 
only possible over ~ 10 MHz in Scheme A. 

To get a better sense of the sacrifices one makes to get a wider bandwidth (due to the diffusion problem) we plot in Fig. llOl a) 
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FIG. 9: Dependence of EIT resonance on gradient and bandwidth. (a,b) The resulting effective susceptibility for Scheme B, with with 
parameters as in Fig. 0(5* = 2 G/mm (black) and S = 8 G/mm (gray)), when one chooses the mismatch such that (Sdisp — 5 , s5 , rC s)wint/2 = 
^(max) -p ne horizontal bars (black and gray) indicate the bandwidth /3 g . (c) The maximum delay r^ max ' = S c s D max /2 for various gradients 
5s versus the resulting bandwidth /3 c ff = SresSsWmt in Scheme A (smaller, blue) and B (larger, red). We have kept the pump power 5 mW 
and the interaction area Wint — 2 cm, hint = 0.5 mm constant, and chosen p — 10 Torr for Scheme A and p = 25 Torr for Scheme B. The 
inset shows a zoom in on the smaller delay points. 
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FIG. 10: Dependence of delay capabilities on field gradient buffer gas pressure, (a) The maximum delay bandwidth product T^ max /3 e g 
versus the slope Sb for the same parameters as in Fig. |9[ c )- The inset shows a zoom in on the lower delay-bandwidth product points. Note 
that Sb =1 G/mm corresponds to /3 e a = (27r)23 MHz. (b) Same plot, but adjusting the width Wi nt so the bandwidth is constant, keeping the 
total pump power at 5 mW. The Scheme B series is chosen so that the total magnetic field difference across Wi nt is 40 G (f3 e g = (27r)56 MHz) 
and for the Scheme A series is 10 G (/3 e ff = (27r)14 MHz), (c) Delay/bandwidth product versus pressure keeping w m t = 2 cm. The solid, 
large (red online) dots are Scheme B for Sb = 2 G/mm and the non-solid (light red) circles are for Sb ~ 4 G/mm. For Scheme A for 
Sb = 0.5 G/mm (small, dark blue dots) and Sb = 1 G/mm (x's, light blue). 



the maximum delay -bandwidth product versus magnetic field gradient (with the same parameters as Fig.|9jc)). We see that for 
Scheme B we maintain full pulse separation capabilities up to bandwidths /3 e g ~ 50 MHz. 

Because of the extremely strong dependence of A gra( j with power we found numerically that one always benefits from us- 
ing smaller h ln t, so this parameter should be chosen to be the smallest reasonable value over which the beam can be easily 
focused. Though too small a value hint would lead to a higher loss from 7diff this problem is almost always dominated by the 
diffusion in the magnetic field gradient and so is not a big consideration. Similarly, in choosing Wi n t we found that the gain in 
transparency from higher intensities tended to outweigh the gain one got from lower slopes. Figure IToF b) shows the maximum 
delay-bandwidth product versus Sb, but keeping the pump power constant and adjusting the interaction with w- m t such that the 
bandwidth /3 e g = (27r)56 MHz was also kept constant. Indeed one sees that one gains by using steeper gradients over smaller 
areas. Ultimately, the slope Sb that can be used in practice will be determined by the manner in which the magnetic field 
gradient and signal dispersion can be generated. 

In the calculations in Figs. I9land ll0f a-b) we chose pressures p = 10 Torr and 25 Torr for Schemes A and B, respectively. 
Numerically these were found to be about optimal. In Fig. [K)] we plot the dependence of the of the maximum delay -bandwidth 
on the pressure for several cases. In Scheme A, the optimal pressures are larger than in the homogenous case due to increased 
importance of reducing diffusion in the magnetic field gradient. Interestingly, even in Scheme B, higher pressures eventually 
reduce the performance. This can be understood from the factor m the analytic estimate for ^4 gra d above. The physical 
origin of this factor is the fact that the EIT width W — H,p/y/8^ e T r decreases with 7 e and therefore makes the resonance more 
sensitive to the averaging over nearby magnetic fields Eq. (I12> . We plot the dependence on pressure for two different gradients 
in each scheme. While the optimal pressures are slightly different, the dependence on pressure is rather weak and so a sensitive 
parameter search versus p should not be required. 



IV. SUMMARY 



We have performed a comprehensive and systematic analysis of EIT resonances, and the resulting pulse propagation char- 
acteristics, in 87 Rb vapors, including effects of couplings to additional levels in the hyperfine structure and a buffer gas. We 
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then calculated the delays, transmissions, and bandwidths for propagation of light tuned to these resonances. We analyzed two 
particular A level schemes (diagrammed in Fig. ^ and found that Scheme B was far superior, in terms of achievable delays 
and delay-bandwidth products, due to the lack of coupling to additional levels. Despite its poorer performance, Scheme A still 
provides reasonable performance and may be desirable since it is much easier to initialize, simply with optical pumping. Im- 
portantly, we found the EIT resonance could be shifted over a wide range of frequencies by applying a homogeneous magnetic 
field, and that the resonance characteristics were quite insensitive to this field over range of about 500 MHz. This analysis serves 
as a useful model to study EIT in conventional slow light, and also as a basis for study of our channelization architecture. 

We then presented a model to analyze the effect of an inhomogenous magnetic field, which causes a strong variation of 
the EIT resonance frequency in the transverse direction. This was then applied to analyze the performance of our proposed 
channelization architecture for wide-band slow light, where a signal pulse is spatially dispersed according to frequency and 
an inhomogeneous magnetic field is applied in a such a way that the EIT resonance frequency matches this dispersion. We 
found that by choosing the magnetic field gradient so the change in the two-photon resonance is slightly mismatched from the 
transverse dispersion of the signal, one could achieve EIT and slow light conditions over a much larger bandwidth than with 
conventional slow light. This is essential for applications in many signal processing applications. We found that the diffusion of 
atoms in the field tended to reduce the delay-bandwidth products with bandwidth. In Scheme B, this architecture should allow a 
delay -bandwidth greater than unity up to bandwidths of about ~ 50 MHz, where delays are ~ 5 ns (see Fig. [8}. Furthermore, 
we note either the pump field power or the magnetic field gradient can be used to control S e s and thus the delay, making it a 
controllable time delay system. 

The buffer gas is important in reducing the diffusion of atoms from into regions of widely varying magnetic field and so higher 
buffer gas pressures are generally desirable for higher magnetic field gradients. However, we also found that higher pressures 
narrow the EIT feature and can therefore increase the sensitivity of the dispersive slope and absorption profile to magnetic field 
gradients. Balancing these two considerations leads to an optimal pressure, which we found this optimum to be near p ~ 10 Torr 
for Scheme A and p ~ 25 Torr for Scheme B, for reasonable parameters. This optimal pressure was not very sensitive to 
the exact value of the gradient and other parameters. We also found that one generally benefited from tight focusing and high 
magnetic gradients. 

In future work, it will be useful to consider the effects of atomic diffusion at a more microscopic level. In particular, it has been 
found that the model used here for diffusion out of the interaction region may overestimate the loss in real systems due to the fact 
that atoms can diffuse back into the interaction region l2lll . Additionally, dynamical jumps of velocity of individual rubidium 
atoms upon collisions with the buffer gas has also been found to be an important consideration fl7ll . Finally, for implementation 
of this system, work is also needed to develop optimal methods for transversely dispersing the signal field and producing large 
linear magnetic field gradients. 

The role played by the differential phase shift in this system is interesting in its own right and merits further investigation. It 
is not entirely clear that the pulse will not be significantly more slowed than our analysis here shows, due to subtleties with the 
transverse dispersion of the signal. Perhaps the signal dispersion or magnetic field could be engineered in such a way that the 
group velocity is governed by the local (and much larger) S, rather than S c g, allowing much larger delays. Furthermore, it may 
be possible to combine this method with aspects of previous light storage experiments I9l [l(jll to significantly increase the delay 
times. 
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